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Abstract. PKS 0537^41, a bright 7-ray emitting blazar, was observed at radio, optical, UV and X-ray frequencies during 
various EGRET pointings, often quasi-simultaneously. In 1995 the object was found in an intense emission state at all wave- 
lengths. BeppoSAX observations made in 1998, non-simultaneously with exposures at other frequencies, allow us to charac- 
terize precisely the spectral shape of the high energy blazar component, which we attribute to inverse Compton scattering. The 
optical-to-7-ray spectral energy distributions at the different epochs show that the 7-ray luminosity dominates the bolometric 
output. This, together with the presence of optical and UV line emission, suggests that, besides the synchrotron self-Compton 
mechanism, the Compton upscattering of photons external to the jet (e.g., in the broad line region) may have a significant role 
for high energy radiation. The multiwavelength variability can be reproduced by changes of the plasma bulk Lorentz factor 
The spectrum secured by lUE in 1995 appears to be partially absorbed shortward of ~1700 A. However, this signature is not 
detected in the HST spectrum taken during a lower state of the source. The presence of intervening absorbers is not supported 
by optical imaging and spectroscopy of the field. 

Key words, galaxies: active — BL Lacertae objects: individual: PKS 0537^41 — ultraviolet: galaxies — X-ray: galaxies — 
Gamma rays: observations 



1. Introduction 

Radiation in the MeV-GeV energy range has been firmly de- 
tected by EGRET in 65 active galactic nuclei, all of blazar 
type (von Montigny et al. 1995; Hartman 1999; Hartman et 
al. 1999). The 7-ray emission, coupled to the small emitting 
volumes inferred by variability time scales, implies beaming in 
a relativistic jet (which makes the emitting regions transparent 
to 7-ray photons, e.g., McBreen 1979; Maraschi et al. 1992), 
a distinctive characteristic of blazars. Among EGRET blazars, 
many exhibit variability on a range of timescales from years to 
days (Wehrle et al. 1998; Mattox et al. 1997; Bloom et al. 1997; 
Hartman 1996; Hartman 1999). 

There is a general consensus that the emission of blazars at 
energies higher than 10 MeV should be attributed to inverse 
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Compton scattering of relativistic electrons off soft photons 
produced in or in the vicinity of the jet. However, the exact role 
played in the scattering process by photons created in the jet 
(synchrotron photons) or outside (broad emission line region 
or accretion disk photons) is not known, nor is the cause of 
the huge, dramatic 7-ray flares observed. Depending on which 
physical parameters are varying (electron density, magnetic 
field density, bulk Lorentz factor), variations of different ampli- 
tude are expected in the synchrotron and inverse Compton com- 
ponents (see Hartman et al. 1996; Wehrle et al. 1998; Ghisellini 
& Madau 1996; GhiselHni & Maraschi 1996). The most ef- 
fective way to discriminate among the different scenarios con- 
sists in observing these sources simultaneously at several wave- 
lengths spanning a broad range. 

PKS 0537-441 (z = 0.896) is one of the most luminous 
and variable blazars at all frequencies, and has been the tar- 
get of monitoring from radio to X-rays at many epochs (Cruz- 
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Gonzalez & Huchra 1984; Maraschi et al. 1985; Tanzi et al. 
1986; O'Brien et al. 1988; Bersanelli et al. 1992; Edelson et 
al. 1992; Falomo et al. 1993a; Treves et al. 1993; Falomo et 
al. 1994; Romero et al. 1994; Heidt & Wagner 1996; Tingay et 
al. 1996; Sefako et al. 2000; Tornikoski et al. 2001). The source 
was detected by EGRET for the first time in 1 99 1 (Michelson et 
al. 1992; Thompson et al. 1993), and then re-observed at many 
successive epochs. It is bright and variable in 7-rays, and its lu- 
minosity in the highest state is comparable to the average lumi- 
nosity of the strongest and best studied EGRET blazar, 3C 279. 
Multiwavelength modeling, based on non-simultaneous 7-ray 
and lower frequencies data, has been proposed by Maraschi et 
al. (1994a) within the synchrolron self-Compton scheme. 

On the basis of its radio variability characteristics (Romero 
et al. 1995) and an off-centered surrounding nebulosity (Stickel 
et al. 1988), it was proposed that PKS 0537^41 is microlensed 
by stars in a foreground galaxy. Search for extended optical 
emission around PKS 0537^41 has a long and controversial 
history (e.g., Falomo et al. 1992; Lewis & Ibata 2000; Scarpa 
et al. 2000). Clarifying the nature of this emission is of impor- 
tance both for the study of the properties of the galaxy and en- 
vironment hosting this very active nucleus and for the alleged 
possibility of microlensing effects. 

In this paper we present multiwavelength observations of 
PKS 0537-441 at various epochs during the EGRET lifetime, 
and particularly focus on the 7-ray flare of 1995 (Sect. 2.1). 
We also report on 1998 BeppoSAX observations in the 0.2-50 
keV band (Sect. 2.2) and consider archival HST spectra (Sect. 
2.3). In Sect. 2.4.2 we compare and discuss the results of all 
imaging studies of the field. In Sect. 3 we discuss the overall 
energy distribution and its implications for the nuclear emission 
mechanisms. 



2. Observations, data analysis and results in 
individual bands 

2. 1. -f-rays 

PKS 0537^41 was repeatedly observed in flie band 0.03-10 
GeV by EGRET aboard CGRO: a summary of the observa- 
tions, with flux levels or upper limits, is reported in the 3rd 
EGRET catalog (Hartman et al. 1999). We have re-analyzed 
the July- August 1991 spectrum with up-to-date software and 
obtain a marginally steeper slope (see Table 1) than found by 
Thompson et al. (1993). During January 1995 the source was 
detected with almost 10-<t significance, with an average flux 
of 90 X 10^*^ photons s^^ cm^^, its brightest recorded 7-ray 
state (see Hartman et al. 1999 for details of the data reduction 
and analysis). The high signal-to-noise ratio of the detection al- 
lowed a rather detailed spectral analysis (see Mukherjee et al. 
1997 and Table 1) and a study of the variability on a day time 
scale. In Fig. 1 we report the EGRET light curve of the blazar 
during the January 1995 pointing, with a temporal binning of 2 
days (see also Hartman 1996). Toward the end of the monitor- 
ing, the flux increases in ~2 days by a factor of at least ~3 with 
respect to the average level during the first 10 days of the obser- 
vation. The target was detected by EGRET also three months 
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Fig. 1. EGRET fight curve of PKS 0537^41 in January 1995 



later, in a state ~4 times lower than shown here (Hartman et al. 
1999). 

In Table 1 and Fig. 2 we report the EGRET data taken 
quasi-simultaneously (within ~6 months) with data at lower 
frequencies. The January 1995 spectrum is significantly flatter 
than measured in 1991, consistent with a common characteris- 
tics of 7-ray blazars, which tend to have flatter spectra during 
flare states (Sreekumar et al. 1996; Sreekumar et al. 2001). 

2.2. X-rays 

X-ray observations of PKS 0537-441 prior to 1991 are summa- 
rized by Treves et al. (1993). The results of die ROSAT 1991 
observation presented by those authors are reported in Table 1 
and Fig. 2. 

We observed PKS 0537-441 with BeppoSAX in 1998 
November 28.3142-30.5291 (UT) as part of a program fo- 
cussed on BL Lacs with flat X-ray spectra, supposedly dom- 
inated by the inverse Compton component. The total on-source 
integration times were 23672, 47396, and 36743 seconds for 
the LEGS (0.1-4 keV), MECS (1.6-10 keV) and PDS (13-300 
keV) instruments, respectively (see Scarsi 1993 and Boella et 
al. 1997 for overviews of the BeppoSAX satelUte). 

The spectra taken by the LEGS and MECS instruments 
have been extracted from the linearized event files using radu of 
8" and 4", respectively, and corrected for background contam- 
ination with library files available at the BeppoSAX Science 
Data Center. The net count rates, averaged over the whole 
pointing, are (1.93±0.11)xl0-2 cts s"^ in the LEGS and 
(2.37±0.08)xlO-2 cts s^^ in the MECS. The spectra mea- 
sured by the PDS instrument have been accumulated from on 
source exposures of the collimator units, and corrected using 
background spectra obtained during the off-source exposures. 
A correction for the energy and temperature dependence of the 
pulse rise time has been also applied (Frontera et al. 1997). The 
average PDS spectrum of the source exhibits significant signal 
up to ~30 keV and appears featureless; its net count rate is 
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Instrument 


a 




V 


Kei. 


1991 Feb 11 


ESO 1.5m-l-CCD 


1.35±0.05 


2.19±0.04mJy 


5.45 X 10" 


1 


1991 Apr 16 


ROSAT-^PSPC 


1.1±0.4 


0.79±0.05 ^Jy 


2.41 X 10^'^ 


2 


1991 Jul 26- Aug 08 


CGRO-l-EGRET 


1.50±0.32 


39±8 pJy 


9.66 X 10^^ 


3,4 


1992 May 14-Jun 04 


CGRO-l-EGRET 




< 48^= pJy 


9.66 X 10^^ 


5 


1992 May 21.41 


IUE-hLWP 




0.66±0.05-'' mJy 


1.15 X 10^^ 


4 


1992 May 21.55 


lUE-^LWP 




0.78±0.03-'^ mJy 


1.15 X 10^^ 


4 


1993 Jul 12 


HST-I-FOS-I-G130H 




0.09±0.01-'^ mJy 


2.14 X 10^^ 


4 


1993 Sep 16 


HST-I-FOS-I-G270H 


1.92 ±0.09 


0.371 ±0.007-'' mJy 


1.15 X 10^^ 


4 


1995 Jan 10-24 


CGRO-l-EGRET 


0.96±0.18 


137±21 pJy 


9.66 X 10^^ 


6 


1995 Jan 30.37 


lUE-l-SWP 


1.2±0.F 


0.90±0.02-'' mJy 


1.67 X 10^^ 


4 


1995 Jan 31.44 


IUE-hLWP 


1.2±0.F 


1.45±0.08-'' mJy 


1.15 X 10^^ 


4 


1995 Feb 01.17 


ESO SEST 




5.16±0.21 Jy 


90 X 10^ 


4 


1995 Feb 03 


ESO 1.5m-l-CCD 


1.29±0.04 


3.03±0.15 mJy 


5.45 X 10" 


7 


1995 Feb 05 


ESO 1.5m-l-CCD 


1.27±0.05 


3.64±0.20 mJy 


5.45 X lO" 


7 


1995 Feb 07 


ESO 1.5m-l-CCD 


1.09±0.04 


5.26±0.25 mJy 


5.45 X 10" 


7 


1995 Feb 27.14 


ESO SEST 




2.88±0.23 Jy 


230 X 10^ 


4 


1998 Nov 28-30 


BeppoSAX 


0.80±0.13 


0.46±0.09 MJy 


2.41 X 10^^ 


4 



" Spectral index {Fv tx /z^""). Uncertainties, both for a and fluxes, are at 90% confidence level 
for the X- and 7-ray measurements, and at 68% for UV, optical and millimetric. No spectral fit has been 
tried when the data signal-to-noise ratio was too low or the spectral range too limited. 

^ Flux density. EGRET flux conversion follows Thompson et al. (1996). Optical-to-X-ray data are corrected 
for Galactic extinction. 

Frequency to which the flux density refers, in Hz. 

References. 1: Falomo et al. 1994; 2: Treves et al. 1993; 3: Thompson et al. 1993; 

4: This paper; 5: Hartman et al. 1999; 6: Mukherjee et al. 1997; 7: Scarpa & Falomo 1997. 

2-(T upper limit. 

^ The uncertainty is only statistical. For lUE data, this was evaluated following Falomo et al. (1993b). 
* From a spectral fit over the band 1700-5500 A. 



(3.9±1.6)xl0 cts s ^ . No emission variability larger than 
15% has been detected within the pointing. 

After rebinning the LECS, MECS and PDS spectra in inter- 
vals where the signal exceeds a 3-(7 significance, we fitted them 
to a single absorbed power-law using the XSPEC routines and 
response files available at the Science Data Center. The LECS 
and PDS spectra normalizations relative to the MECS have 
been treated as a free (best fit value 0.86±0.09, within the ex- 
pected range) and fixed (0.85) parameters, respectively, and the 
neutral hydrogen column density has been fixed to its Galactic 
value (A'^Hi = 2.91 x lO^" cm^^ Murphy et al. 1996). The fit 
is satisfactory (x^ — 0.6). The best-fit power-law parameters 
are reported in Table 1 and the deconvolved de-absorbed X-ray 
fluxes in the BeppoSAX energy band are reported in Fig. 2. We 
also searched for a possible iron emission line coiTesponding 
to the Ka transition at the expected redshifted energy on the 
MECS spectrum, but we do not detect any line with intensity 
larger than 1.4 x 10^^'* erg s^^ cm^^ (3ct upper limit). 

The spectrum and emission state detected by BeppoSAX 
are consistent with those measured at earlier epochs by 
EXOSAT and Einstein (Treves et al. 1993), while the 
BeppoSAX flux at 1 keV is almost a factor of ^2 lower than 
that detected by ROSAT in 1991. 



2.3. Ultraviolet 

lUE performed low dispersion spectroscopy of PKS 0537- 
441 in the 2000-3000 A (LWR or LWP cameras) and 1200- 
1950 A (SWP) wavelength ranges at various epochs between 
1980 and 1995, and in two occasions. May 1992 and January 
1995, quasi-simultaneously with EGRET pointings (Table 1). 
The coiTesponding lUE spectra were retrieved from the INES[] 
and NEWSIPS| (Nichols & Linsky 1996; Garhart et al. 1997) 
archives and, for the 1995 data only, extracted from the bi- 
dimensional spectral images with the GEX routine (Urry & 
Reichert 1988). 

The INES and NEWSIPS spectra of May 1992 agree in flux 
to within ^10%, and therefore can be considered consistent, 
given the low level of the source and the possible differences 
in background estimate between the two extraction methods. In 
Table 1 and Fig. 2 we have reported the results obtained with 
the NEWSIPS extraction. As a consistency check, we note that 
these UV fluxes are compatible with the power-law fitting the 
optical data taken more than one year earlier (see Table 1 and 
Fig. 2). 

For the spectra of January 1995, the INES and NEWSIPS 
methods return compatible results for the lUE SWP range (the 

' http://ines.vilspa.esa.es/ines/docs/contents.html 
^ http://archive.stsci.edu/iue/index.html 
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Fig. 2. Radio-to-7-ray spectral energy distribution of PKS 0537^41 at different states (see Table 1 for data references). The filled 

pentagons, triangles and circle represent millimetric, optical-to-UV and 7-ray data in Jan-Feb 1995, respectively; diamonds are 
IRAS 1984 data from Impey & Neugebauer (1988). The open triangles stand for optical February 1991 and UV 1992 data, open 
squares represent ROSAT April 1991 (large symbol) and BeppoSAX 1998 (small symbol) data, and open circles represent the 
EGRET Jul- Aug 1991 detection and 1992 upper limit. Data at optical-to-X-ray frequencies are corrected for Galactic extinction 
(see text). Superimposed on the multiwavelength data are the model curves for the low 1991-1992 state (dotted), and for the high 
1995 state (solid). The model parameters for the low state are 7min = 1, 7b = 400, 7max = 2.5 x 10^, ni = 1.6, n2 = 4, 
B = 3.5 Gauss, K = 6 x 10^ cm^^, R = 2.7 x lO^^ cm, 5 ~ F = 10 (see Ballo et al. 2002 for the parameter notation). In 
high state, the above parameters remain unchanged, except that F = 11 and n2 = 3.8. The external Compton seed photon source 
(disk photons reprocessed in the broad line region) is assumed to have a density ?7ext = 6 x 10~^ erg cm~^; the disk temperature 
is assumed to be 10^ K 



flux at 1800 A, averaged between the two extractions and cor- 
rected for Galactic extinction, is 1.21 ±0.03 mJy), while in the 
LWP range the NEWSIPS extracted flux is up to a factor of 2 
higher than found wifli INES. The GEX extracted flux is ~25% 
lower than that given by INES and NEWSIPS in the SWP band, 
and is consistent with the INES result in the LWP band. Since 
we have verified that the differences between pairs of lUE 
spectra of bright sources (e.g., the BL Lac object PKS 2155- 
304) extracted with different methods do not exceed 1-2% (see 
also Urry et al. 1993), we ascribe the discrepancies we find 
for PKS 0537-441 to the difficulty of correctly evaluating the 
background at UV flux levels comparable with the lUE sensi- 
tivity limit and of properly removing the solar scattered light, 
which critically contaminates the lUE LWP camera acquisi- 
tions after 1992 (CapUnger 1995, and references therein). 



UnUke tiie NEWSIPS and INES spectra, the fluxes result- 
ing from the GEX spectra of PKS 0537^41 are consistent 
with the power-law which fits the quasi-simultaneous (within 
few days) optical data (see also below). Therefore, we have de- 
cided to adopt the GEX extraction for the 1995 lUE spectra. 
These have been cahbrated according to Bohhn et al. (1990, 
SWP) and Cassatella et al. (1992, LWP), and have been binned 
in 70-100 A-wide wavelength intervals, after removing cos- 
mic rays and by avoiding regions affected by spurious sunUght 
(A > 2800 A) and by large noise (1900 A < A < 2300 A). The 
binned spectra are reported in Fig. 3. 

No emission lines superimposed on the UV continuum are 
detected in any of the lUE spectra. In particular, the Lya and 
C IV lines, well detected by HST (see below), fall at the edges 
of the LWP range, where the sensitivity is poor. Moreover, the 
high state detected by lUE in 1995 (a factor of ~4 higher than 
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Fig. 3. lUE (circles) and HST FOS (solid curves) spectra of 
PKS 0537-441 not corrected for Galactic extinction. The lUE 
data refer to a bright (January 1995, filled circles) and a faint 
(April 1982, open circles) state. The FOS G130H and G270H 
spectra, smoothed over wavelength intervals of 12.5 and 5 A, 
respectively, sample a low state (July and September 1993, re- 
spectively), comparable to that observed by lUE in 1982. Note 
that the strongest emission lines, Lya and C IV, well visible in 
the HST spectrum, occur out of the usable wavelength region 
of the lUEdata 



observed by HST), may have implied smaller line equivalent 
widths. A suggestion of C IV emission may be present in the 
average of the LWP spectra of PKS 0537^41 taken before 
1988 (Courvoisier & Paltani 1992). 

A power-law with slope Oiv ^ 1.2 (/^ cx v^""") fits 
marginally well (reduced — 2) the 1995 optical and UV data 
at wavelengths longer than ~1700 A, corrected for Galactic 
dust extinction (£^b-v = 0.037 mag, Schlegel et al. 1998) with 
the law by Cardelli et al. (1989). The spectral index is in agree- 
ment with that determined for the optical spectra only (Scarpa 
& Falomo 1997). The LWP spectrum appears steeper than this, 
which may account for the not completely satisfactory value of 
the . The discrepancy may be due to the limited quality of 
the lUE LWP data, and to the non-strict simultaneity. 

The UV points at wavelengths shorter than 1700 A clearly 
deviate from the optical-to-UV power-law behavior and sug- 
gest a flux deficit by up to 50% at -1350 A (Fig. 3). This fea- 
ture is seen independent of the method adopted for spectral ex- 
traction. However, we note that the background has an irregular 
shape along the dispersion direction, and possible instrumen- 
tal effects cannot be excluded. The only other SWP spectrum 
available for PKS 0537-441 (25 April 1982, Fig. 3) obtained at 
a very faint state, has insufficient signal-to-noise ratio to con- 
firm or disprove the reality of the feature (see also Lanzetta et 
al. 1995). Moreover, a similar dip is not seen in the HST FOS 
G130H spectrum (see below), which has a level comparable 
to that of the 1982 lUE spectrum (Fig. 3), but a much better 
signal-to-noise ratio. Therefore, we do not consider this fea- 



Fig.4. FOS spectrum (G270H) of the blazar PKS 0537^41. 
Strong nuclear emission lines (z = 0.896) and faint Galactic 
absorption features are clearly detected 



ture in the analysis of the multiwavelength energy distribution 
(Sect. 3). 

PKS 0537^41 was observed by the HST Faint Object 
Spectrograph (FOS) in July and September 1993 (pre- 
COSTAR) with the G130H and G270H gratings, respectively 
(Table 1). We retrieved both spectra from the HST archive| 
(see also Bechtold et al. 2002). The G130H spectrum, which is 
rather noisy, presents a strong, sharp emission feature at ^^1300 
A that cannot be identified with any known emission (if at rest 
frame A — 690 A). Because of this and the unusual shape 
we believe it is an instrumental effect. On the other hand, the 
G270H spectrum has a good signal-to-noise ratio, and exhibits 
broad, intense emission lines and Galactic interstellar absorp- 
tions (see Fig. 4). 

We derived the central wavelength of each emission line by 
fitting the core of the features with Gaussians. Since the Lya 
profile is slightly asymmetric due to a possible narrow absorp- 
tion on the blue side and close to the center of the line, we 
derived the centroid of this line by excluding the inner por- 
tion of the core. On the basis of their measured central wave- 
lengths (see Table 2), the emission lines are consistently iden- 
tified with Lya, Si IV, and C IV at z = 0.896 ± 0.001. This 
is slightly larger than, but still consistent with historical red- 
shift determination (z ~ 0.894, Peterson et al. 1976), based 
on Mg IIA2798 emission, but only marginally consistent with 
z = 0.892 ± 0.001 given by Lewis and Ibata (2000), based 
on the Balmer line series and O III. The intensities and equiv- 
alent widths of the UV emission lines are reported in Table 2. 
After excluding the spectral regions affected by the lines, the 
UV continuum is well fitted by a power-law of index a ~ 1.9. 

We computed the continuum fluxes from the HST and lUE 
spectra by averaging the signal in 200-A-wide bands in spec- 
tral regions free from lines, spurious features and noise, and 
corrected for Galactic extinction (see Table 1 and Fig. 2). The 



http://archive.eso.org/wdb/wdb/hst/science/form 
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Table 2. UV Emission Lines of PKS 0537^41 



Ion 


A" 


EW" 


r 


Lya A1216 


2305.9 ± 0.5 


30 ±5 


5.5 ±0.8 


Si IV A1403 


2660.0 ± 0.5 


1.5 ±0.2 


0.25 ± 0.04 


CIVA1549 


2937.0 ± 0.5 


22 ±3 


3.6 ±0.5 



" Observed line wavelength in A. 
Observed equivalent width in A. 
Line intensity in 10~^* erg s~^ cm~^, corrected for 
Galactic extinction. 



UV state of 1995 was about a factor 2 higher than in 1992 and 
a factor 4 higher than in 1993. The optical-to-UV spectrum in 
1995 was flatter than the HST UV spectrum in 1993, consistent 
with the expected correlation of spectral flattening with bright- 
ening in synchrotron sources. 

2.4. Optical 

2.4.1 . Spectroscopy 

Optical CCD spectrophotometry was obtained in February 
1991 (Falomo et al. 1994) and on 3,5 and 7 February 1995 at 
the 1.5m telescope of the European Southern Observatory (see 
Scarpa & Falomo 1997 for details of data acquisition, reduction 
and analysis, and Table 1 and Fig. 2 for the results). A 70% in- 
crease is apparent between 1991 and 1995 in the optical flux. 
We detected a monotonic increase of about 70% during the 3 
nights of the 1995 monitoring, accompanied by spectral hard- 
ening. The Mg IIA2798 emission line is observed in all spec- 
tra, with non- significantly variable intensity, Ix c::^ 6 x 10~^^ 
erg s~^ cm~^ (not corrected for Galactic reddening, Scarpa & 
Falomo 1997; Treves et al. 1993. Note that the line intensity in 
the latter reference is misreported by a factor of 10). 

2.4.2. Imaging 

Stickel et al. (1988) reported the detection of extended opti- 
cal emission surrounding PKS 0537^41, the characteristics of 
which led them to interpret it as a foreground disk galaxy along 
the line of sight to the blazar. However, this finding was not 
confirmed either by higher resolution optical images (Falomo 
et al. 1992) and near-IR imaging (Kotilainen et al. 1998), con- 
sistent with the high redshift of the source and the relatively 
bright nucleus. This source was also imaged by HST using 
WFPC2 and F702W filter (Urry et al. 2000; Scarpa et al. 2000) 
and found to be unresolved. However, based on the analysis of 
the same HST images, Lewis and Ibata (2000) claim the find- 
ing of faint non-axisymmetric extended emission. This result 
was obtained from the subtraction of a scaled PSF constructed 
with the Tiny Tim software (Krist 1995). While Tiny Tim offers 
a reliable description of the central part of the HST-WFPC2 
PSF it does not account for the emission due to the large an- 
gle scattered Ught (particularly for the Planetary Camera) that 
produces extra emission at radii larger than ~2 arcsec, as il- 
lustrated in Fig. 5 (see also Scarpa et al. 2000). We therefore 
beUeve that the extra Ught reported by Lewis and Ibata (2000) 



is fully accounted for by the contribution of the scattered Ught. 
Moreover, we note that no signature of stellar features has been 
reported in the optical spectra of PKS 0537-441 to support the 
presence of a foreground galaxy (Stickel et al. 1988; Lewis & 
Ibata 2000). Therefore, the microlensing hypothesis does not 
appear tenable (see also Heidt et al. 2002). 

2.5. Millimeter 

PKS 0537^41 was monitored during flie year 1995 wifli 
the 15-metre Swedish-ESO SubmilUmetre Telescope (SEST) 
at La Silla, Chile (Tomikoski et al. 2002, in preparation). 
Observations were made at 90 GHz using a dual polariza- 
tion Schottky receiver, and at 230 GHz using a single channel 
bolometer. More details on the observations and data reduc- 
tion can be found in Tornikoski et al. (1996). The sampling 
is only marginally simultaneous with EGRET. The fluxes ob- 
tained closest in time with the EGRET observation are reported 
in Table 1, along with errors, which consist of the standard de- 
viation among the individual daily measurements and the esti- 
mated effect of the caUbration and pointing uncertainties. 

These data suggest a very moderate brightening in the radio 
flux curve at 90 GHz just about the time of the multiwavelength 
campaign, but not comparable to the 1994 outburst (Tornikoski 
et al. 1996). In the 230 GHz data no changes can be seen from 
December 1994 to the end of February 1995, but the sampUng 
is too sparse to completely exclude the possibility of a minor 
burst similar to that seen in the 90 GHz curve (see Tornikoski 
et al. 1996). 

3. Discussion 

In Fig. 2 we have reported the broad-band spectral en- 
ergy distributions of PKS 0537^41 constructed with quasi- 
simultaneous data from the millimetric to the 7-ray frequen- 
cies (see Table 1). In addition, we have reported the un- 
pubUshed BeppoSAX spectrum. The multiwavelength state 
of PKS 0537^41 in January-February 1995 was one of the 
brightest recorded for this source during the lifetime of EGRET 
(cf. Tornikoski et al. 1996; Falomo et al. 1994; Hartman et al. 
1999). The comparison of 7-ray detection and optical flux of 
1991 and the EGRET upper limit and UV flux of 1992 is sug- 
gestive of Uttle or no multiwavelength variability between the 
two epochs. The increase of the optical-to-UV flux from the 
1991-1992 to the 1995 level corresponds to a variation of sim- 
ilar amplitude (a factor of ~2) in the 7-rays, or only sUghtly 
larger. The flat BeppoSAX spectrum suggests that a single 
emission component dominates in the energy band 0. 1-30 keV. 

PKS 0537^41 exhibits the typical double-humped multi- 
wavelength spectral shape of "Low-frequency-peaked" blazars 
(Padovani & Giommi 1995; Sambruna et al. 1996; Fossati et 
al. 1998): the first component peaks at wavelengths longer than 
the optical (likely in the far-infrared, as suggested by the IRAS 
data, taken at a much earlier epoch, Impey & Neugebauer 1988) 
and is due to synchrotron radiation. The second component, 
which peaks around ~ 10^^ — 10^^ Hz (Fig. 2) and is a factor 
^6 more powerful than the synchrotron maximum, is proba- 
bly produced via inverse Compton scattering between relativis- 
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Fig. 5. The upper panels show, from left to right, the HST Planetary Camera images in F702W filter of PKS 0537-441, the 
unresolved blazar H 1722+1 19, and a star (the angular size of each field is 9". 2 x 9". 2). The lower panels report, from left to 
right, the PSF model fit (soUd curve), including the contribution of scattered light, to the observed profiles of surface brightness 
(in mag arcsec"^) vs radius (in arcsec) for PKS 0537-441, H 1722+1 19, and many stars (dots). In the third lower panel we also 
report, for comparison, the Tiny Tim PSF model (dashed curve), which is clearly unable to describe the profile of point-Uke 
sources at radii larger than ~2 arcsec (see also Scarpa et al. 2000) 



tic electrons and synchrotron photons (Maraschi et al. 1992; 
Maraschi et al. 1994b; Bloom & Marscher 1993) or external 
photons (broad emission line region or accretion disk, Dermer 
& Schlickeiser 1993; Sikora et al. 1994). 

The overall energy distribution of PKS 0537^41, the 
inverse-Compton dominance, and the presence of optical and 
UV emission lines suggest that external Compton upscat- 
tering may be significant with respect to the synchrotron 
self-Compton process in producing the high energy emission 
(Ghisellini et al. 1998; Ghisellini 2001). To model the broad- 
band energy distribution of PKS 0537^41 we have assumed 
that the emission is produced in a region filled by relativis- 
tic particles which radiate at low energies via synchrotron, and 
upscatter both synchrotron photons and accretion disk photons 
reprocessed in the broad emission line region. We have not in- 
cluded the direct contribution of the accretion disk to the popu- 
lation of seed photons to be upscattered by the jet electrons, be- 
cause those would be highly redshifted in the frame of the blob, 
which is emitting at a distance of ^ 0.3 pc (see caption to Fig. 
2) from the jet base (Sikora et al. 1994). The size of the emit- 
ting region has been constrained with the variability time scale 
in 7-rays (^^2 days); the electron energy distribution is mod- 
eled with a double power-law (see Tavecchio et al. 2000 and 
Ballo et al. 2002 for more details of the model). The low state 
model has been constrained with the BeppoSAX spectrum and 
with the 1991-1992 optical and 7-ray data, while the high state 
model reproduces the 1995 multiwavelength data. Note that the 
ROSAT spectral measurement of April 1991, which is affected 
by a large uncertainty, is consistent both with the low and high 
state model curves. The IRAS data are marginally consistent 
with the high state model. 

The difference between the two model curves (reported 
in Fig. 2; see caption for the model parameters) is solely de- 
termined by a change in the Lorentz factor of the relativistic 
plasma bulk motion (increasing by 10% from the low to the 
high state) and in the index of the upper branch of the elec- 
tron energy distribution (n2 in the notation of Ballo et al. 2002, 
sUghtly flatter in the brighter state, see caption to Fig. 2). Under 



our assumption that the high energy component results from 
the contribution of both synchrotron- self Compton and exter- 
nal Compton mechanisms, the variation of the bulk Lorentz 
factor would produce a correlation between the variability am- 
plitude of the synchrotron and inverse Compton components 
which is intermediate between linear and half-cubic (Ghisellini 
& Maraschi 1996). This is consistent with the observations, al- 
though the difficulty of exactly locating the peaks of the emis- 
sion components, the non-strict simultaneity of the optical-to- 
UV spectra and the rapid 7-ray flare seen at the end of the 
EGRET viewing period make the fit results only indicative. 

From our measured UV line intensities reported in Table 2 
and from dereddenedMgll emission fine intensity (Sect. 2.4. 1), 
assuming Hq = 65 km s~^ Mpc~^, fim = 0.3, Qa = 0.7, we 
derive a total line luminosity Lblr ~ 4.7 x 10^'' erg s~^ (the 
other observed optical emission lines are sufficiently weak that 
their contribution to the line luminosity is not significant, see 
Lewis & Ibata 2000). Using this and the external photon den- 
sity assumed in our model, J7ext = 6 x 10~^ erg cm~^, we can 
evaluate the size of the broad line region, Rblr — 7.9 x 10^^ 
cm. As an independent check, we have considered the empiri- 
cal relationship determined by Kaspi et al. (2000) between the 
size of the broad line region and the luminosity of the ther- 
mal continuum at 5100 A in quasars. Assuming a covering fac- 
tor of the broad Une clouds of 10%, we can estimate the disk 
luminosity to be idisk ~ 4.7 x 10^^ erg s~^. This value is 
consistent with the approximate upper limit which can be de- 
rived from the broad band spectrum, ~ 2 x 10^^ erg s~^. After 
accounting for a factor of ~3 difference between a bolomet- 
ric and a monochromatic disk output, we find that our disk lu- 
minosity estimate implies, according to Kaspi et al.'s formula, 
-Rblr — 5.8 x 10^'' cm, well consistent, given the uncertain- 
ties, with the size we have determined based on the observed 
I/BLR and on the assumed photon density. We finally note that 
our observed Lblr is about a factor 4 lower than that esti- 
mated for this source by Celotti et al. (1997), based on an old 
measurement of the Mgll emission only. 
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The validity and origin of the spectral dip seen in the 
lUE spectrum of January 1995 at wavelengths shorter than 
1700 A remain to be established. Assuming it to be real, 
it could be qualitatively described by an absorption edge or 
broad trough, and may be identified with neutral hydrogen 
ionization discontinuity (Lyman limit) at the redshift of the 
source or with Lya forest blanketing up to a maximum red- 
shift of ^0.4. Continuum decrements shortward of Lya emis- 
sion, Lyman limit systems and damped Lya absorption features 
have been reported in few blazars at redshifts 0.5 ^ z ^ 1.5 
(e.g., PKS 0637-752, Cristiani et al. 1993; PKS 0735H-178, 
PKS 2223-052 Courvoisier & Paltani 1992, Lanzetta et al. 
1995). Absorption features have recently been reported also in 
X-ray blazar spectra (Tavecchio et al. 2000). The host galaxy or 
the halos of the galaxies located in the vicinity of PKS 0537- 
441 and at similar redshift (projected distances of ~30 kpc, 
Heidt et al. 2002) may be responsible for absorption in the 
bluer part of the UV nuclear spectrum. Using the photoelec- 
tric cross-section computed by Rumph et al. (1994) to model 
the interstellar opacity at extreme UV wavelengths, we esti- 
mated that in the case of Lyman continuum absorption the 
equivalent intrinsic A'hi would be about 3 orders of magni- 
tude lower than the Galactic A'hi in the direction of the blazar, 
and therefore its effect on the X-ray spectrum would be un- 
detectable. Alternatively, the FOS G130H spectrum may sug- 
gest resemblance with those of Broad Absorption Line quasars 
(e.g., Hamann & Ferland 1999; Arav et al. 2001), so that, in 
the lower resolution lUE SWP spectrum, many broad features 
would appear blended in a unique edge. 

Higher signal-to-noise ratio data (ideally HST STIS spectra 
simultaneously covering the wavelength range 1100-2700 A) 
would be necessary to confirm the presence of the feature and 
accurately model it. The nuclear emission and environmental 
characteristics of the blazar PKS 0537^41 make it a case study 
for future 7-ray missions hke AGILE and GLAST and for the 
NGST, respectively. 
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